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ABSTRACT 

X-ray surface brightness fluctuations in the core of the Perseus Cluster are analyzed, 
using deep observations with the Chandra observatory. The amplitude of gas density 
fluctuations on different scales is measured in a set of radial annuli. It varies from 8 
to 12 per cent on scales of ~ 10 — 30 kpc within radii of 30 — 160 kpc from the cluster 
center and from 9 to 7 per cent on scales of ^ 20 — 30 kpc in an outer, 160 — 220 kpc 
annulus. Using a statistical linear relation between the observed amplitude of density 
fluctuations and predicted velocity, the characteristic velocity of gas motions on each 
scale is calculated. The typical amplitudes of the velocity outside the central 30 kpc 
region are 90 — 140 km/s on ^ 20 — 30 kpc scales and 70 — 100 km/s on smaller scales 
^7—10 kpc. The velocity power spectrum is consistent with cascade of turbulence and 
its slope is in a broad agreement with the slope for canonical Kolmogorov turbulence. 
The gas clumping factor estimated from the power spectrum of the density fluctuations 
is lower than 7 — 8 per cent for radii ^ 30 — 220 kpc from the center, leading to a 
density bias of less than 3 — 4 per cent in the cluster core. Uncertainties of the analysis 
are examined and discussed. Future measurements of the gas velocities with the Astro- 
H , Athena and Smart-X observatories will directly measure the gas density-velocity 
perturbation relation and further reduce systematic uncertainties in these quantities. 

Key words: turbulence-methods: observational-methods: statistical-techniques: im¬ 
age processing-galaxies: clusters: intracluster medium-X-rays: galaxies: clusters 


1 INTRODUCTION 

The Perseus Cluster is the brightest galaxy cluster in the 
X-ray sky and has been attracting astronomers attention 


* zhur@stanford.edu 


over the last 40 years. Being the X-ray brightest cluster, 
Perseus has led to many fundamental cluster discoveries, 
such as peaked cluster emission (Fabian et ah; 1974), strong 
FeXXV and FeXXVI lines and thermal emission mechanism 
of the gas ( Mitchell et al.||1976 ), presence of bubbles of rel¬ 
ativistic plasma (Boehringer et al. 1993). The cluster was 
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recently observed deeply with the Chandra (see e.g. Fabian 

et al.|2000, 2011; Schmidt, Fabian, & Sanders|2002; Sanders 

V Fabian||2007 and references therein), XMM-Newton (see 

e.g. Bohringer et al.|2002; |Churazov et al. 2003, 2004; Mat- 

sushita et al.||2013| and references therein) and Suzaku (see 

e.g. Simionescu et al.||2011| Ueda et al. |2013| Werner et al.] 

2013| Tamura et al.||2014| Urban et al. 2014, and references 

therein) observatories, as well as by the previous-generation 

X-ray observatories (see e.g.|Forman et al. 1972| Mushotzky 

et al. 1981; Branduardi-Raymont et al. |1981; Arnaud et al. 

1994| Fabian et al. 1994; |Furusho et al. 2001; Gastaldello 

V Molendi| 2004|. These observations unveiled a variety of 

structures present in the hot gas of the cluster, reflecting 
richness and complexity of different physical processes oc¬ 
curring in the intracluster plasma. For example, an east-west 

asymmetry in the X-ray surface brightness (SB) (|Schwarz 

et al. 1992), aligned with a chain of bright galaxies, sug- 
gests an ongoing modest merger, which can induce sloshing 

of the gas and drive gas turbulence (Churazov et al. |2003|, 

while the central few arcmin are dominated by AGN activ¬ 
ity in a form of inflated bubbles of relativistic plasma (see 

e.g. Boehringer et al. 1993| |Churazov et al.|2000; Fabian et| 

al. 2000) and weak shocks around them. Rising buoyantly 
and expanding, the bubbles uplift cooler X-ray gas from the 
core, producing filamentary structures seen in the optical, 
far-ultraviolet and soft X-ray, and excite internal waves, en- 


ergy of which is transported to gas turbulence (Churazov et 


al.l 

2001 

Omma et al. 2004; 

Forman et al. 2007; 

Fabian et 

al. 

2008 

(Canning et al.||2014 

Hillel & Sokei'|2014 

)• 


The truly spectacular statistics accumulated by the 
Chandra observations of Perseus (over 85 million counts 
within 6 arcmin from the center, with 100s and 1000s counts 
per sq. arcsec) makes this data set an extremely powerful 
tool to probe structures in the intracluster medium (ICM) 
on a range of spatial scales. Carefully modeling the Pois¬ 
son noise, one can probe small-scale fluctuations, which we 
cannot see in the cluster images by visual inspection. Here, 
we present statistical analysis of these SB and density fluc¬ 
tuations, using the power spectrum (PS) statistics. We will 
examine the radial variations of the PS within the Perseus 
core (220 x 220 kpc), while its energy dependence will be 
considered in our future works. 


Statistical analysis of the X-ray SB and/or density, pres¬ 
sure fluctuations was done for two clusters so far: the Coma 
Cluster (Schuecker et al. 2004 Churazov et al. 2012) and 
AWM7 ( Sanders &; Fabian|2012 ). It was shown, for example, 
that the amplitude of density fluctuations in Coma ranges 
from 5 to 10 per cent on scales 30 — 500 kpc, leading to non¬ 
trivial constraints on perturbations of the cluster gravita¬ 
tional potential, turbulence, entropy variations and metal- 
licity. The fluctuations appear to be correlated, implying 
suppression of strong isotropic turbulence and conduction 
( Sanders et al.||2013 ). 


The Chandra data on Perseus allow us to probe fluctua¬ 
tions on similar or even smaller scales, down to few kpc. The 
main questions we address here include: what is the ampli¬ 
tude of density fluctuations and how does it vary with the 
distance from the cluster center (Section |5.1| ); what is the 
power spectrum of gas turbulent motions (Section |5.2| ; how 


clumpy is the gas in the core (Section 5.3). To answer these 
questions, various instrumental effects (e.g. PSF variations, 
detector QE fluctuations), contribution of the unresolved 
point sources, subtraction of the Poisson noise as well as 
its fluctuations need to be treated carefully. We extensively 
investigate various systematic uncertainties of the analysis 
and discuss them in Section [6] 

This paper is complementing our series of publications 
on SB and density fluctuations in the core of the Perseus 
Cluster. One application of the velocities of gas motions 
measured here is discussed in Zhuravleva et al. (2014b). 
There we show that turbulent dissipation in the cluster core 
produces enough heat to balance locally radiative energy 
losses from the ICM at each radius. Therefore, turbulence 
may play an important role in AGN feebdack loop and be 
the key element in resolving the gas cooling flow problem 
( Zhuravleva et al.|2014b ). In our next paper, we will discuss 
the nature of the fluctuations by measuring cross-spectra 
of fluctuations in different energy bands (Zhuravleva et al., 
2015b, in prep.). Similar analysis for the Virgo and Coma 
clusters are in progress. 

Throughout the paper we adopt ACDM cosmology with 
Qm — 0.3, Qa = 0.7 and h — 0.72. At the redshift of 
Perseus, z = 0.01755, the corresponding angular diameter 
distance is 71.5 Mpc and 1 arcmin corresponds to a physical 
scale 20.81 kpc. Position of the cluster center is taken from 
the NASA/IPAC Extragalactic Databas^Jand corresponds 
to the position of the central galaxy NGC1275. The total 
galactic HI column densit^ is 1.36T0 21 cm -2 . Solar abun¬ 
dances of heavy elements are taken from|Anders & Grevesse 
( 1989| ). We define the characteristic scale-dependent ampli¬ 
tude A 3 o(k) through the power spectrum P^Dik) as follows 
A\ D (k) — f 47t P 3 D(k)k 2 dk = 47vP3D(k)k 3 . A wavenumber k 
is related to a spatial scale l without a factor 27r, i.e. k — 1/1. 


2 INITIAL DATA PROCESSING 


For our analysis we use public Chandra data with the fol¬ 
lowing ObsIDs: 3209, 4289, 4946, 4947, 4948, 4949, 4950, 
4951, 4952, 4953, 6139, 6145, 6146, 11713, 11714, 11715, 
11716, 12025, 12033, 12036, 12037. The initial data process¬ 
ing is done using the latest calibration data and following 


the standard procedure described in Vikhlinin et al. (2005). 


For the analysis, we choose 0.5 — 3.5 keV energy band. In 
this band, the gas emissivity has a weak temperature de¬ 
pendence as shown in Fig. [l] and, therefore, the uncertainty 
associated with the conversion of PS of SB fluctuations to 


the PS of density fluctuations is reduced (see Section 5.1). 

Correcting for the exposure, vignetting effect and sub¬ 
tracting the background, a mosaic image of the Perseus Clus- 


1 https://ned.ipac.caltech.edu/ 

2 We use the mean value of the Galactic absorption within the 
field of the Perseus Cluster. Exclusion of soft photons (0.5 — 1 keV) 
from the SB analysis does not change the resulting amplitude of 
the SB and density fluctuations. Therefore, even if fluctuations of 
the Galactic absorption are present, their role in the cluster core 
is minor. 
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Figure 1 . Emissivity of an optically thin plasma as a function 
of temperature observed by the ACIS-S detector on Chandra in 
two energy bands: 0.5 — 3.5 (solid) and 3.5 — 7.5 (dashed) keV. 
Abundance of heavy elements is assumed 0.5 relative to Solar. 
Shaded region shows a range of gas temperatures in the core of 
the Perseus Cluster (3 — 6.5 keV). Notice, that for the 0.5 —3.5 keV 
band, the emissivity is almost independent of the temperature 
within this region. Therefore, the X-ray surface brightness in this 
shaded band lx is just a function of gas electron number density 
n e , namely, lx ~ f ri*dl, where l is the length of the line of sight. 


ter is produced, which is shown in Fig. [2] The combined ex¬ 
posure map is shown in Fig. [3] The total exposure of the 
cleaned data set is ~ 1.4 Ms. Notice, that the exposure cov¬ 
erage is inhomogeneous. It varies from few-10 6 s in the center 
down to few-10 5 s in some distal regions. The uncertainties 
associated with such inhomogeneous coverage are discussed 
in Section [6] 

The Perseus dataset includes observations with differ¬ 
ent offsets. For the analysis of SB fluctuations, it is essential 
to know the PSF at each point within the cluster image. 
We randomly generate positions, where the PSF is checked. 
The PSF at each position for each observation is determined 
from the Chandra PSF libraries (Karovska et al. 2001). By 
combining the PSF model images with weights proportional 
to exposure of each observation, the final PSF map is ob¬ 
tained. These maps are used to correct the PS of the SB 
fluctuations (Section [4). In order to illustrate the variations 
of the PSF shape within the cluster mosaic image, we show 
the PSF map divided by the exposure map in Fig. [3] No¬ 
tice a non-monotonic behavior of the PSF and its elongated 
(distorted) shape close to the edges of the image. 

As the next step, we excise point sources from the im¬ 
age of the cluster. Using the wvdecomp tool, the point source 
candidates are found as peaks in the image with S/N > 4. In 
order to verify the significance of each source detection, we 
smooth heavily the image of the cluster with a 20 /x Gaussian, 
excluding the brightest point sources first. We then subtract 


the smoothed image from the initial image to remove the 
cluster contribution and calculate the point source fluxes in 
small regions around each point source candidate account¬ 
ing for the combined PSF as F p . a . = ^ I p . s .PSF/ ^ PSF 2 , 
where / p . s . refers to the residual image. The uncertainties on 
the fluxes are <rp p s = CiPSF 2 /^ PSF 2 , where Ci is 
the number of counts in the initial cluster image. Finally, 
sources with F p , s Jcjf v . s , >4 are cut out from the images, 
using circles with the radius 1.5 times the 90 per cent radius 
of the combined PSF (see Fig.[2|. Contribution of the unre¬ 
solved point sources is estimated using sensitivity maps and 
Log N - Log S distribution of resolved sources. In case of the 
Perseus core, the contribution of unresolved point sources is 
negligible. 


3 SPHERICALLY-SYMMETRIC MODEL 


The spherically-symmetric radial surface brightness profile 
of the Perseus is shown in Fig. [4] Fitting this profile with 
a model assuming photon counting noise only, gives the 
core radius r c = 1.26 arcmin and the slope of the model 
1 3 = 0.53 with negligible uncertainties. Allowing for 10 per 
cent systematic error on the SB in each annulus, the un¬ 
certainties become 0.1 and 0.01 for r c and fd respectively. 
Throughout the paper, this model is used as a default model 
of the unperturbed cluster. Sensitivity of our results to the 
choice of the model are discussed in Section |6] 

The deprojected thermodynamic properties of the clus¬ 
ter, namely, the number electron density n e and the electron 
temperature T e , are obtained from projected spectra using 
a procedure described in Churazov et al. (2003). The de¬ 


projected spectra are fitted in the broad 0.6 — 9 keV band 


using the XSPEC (Smith et al. 2001 Foster et al. 2012) 


code and APEC plasma model based on ATOMDB version 
2.0.1. The abundance of heavy elements is assumed to be 
either a) constant with radius 0.5 relative to Solar, or b) a 
free parameter in the plasma model. Fig. [4] shows the de¬ 
projected data and their approximations with continuous, 
smooth functions. Using these approximations, the sound 
speed c s for an ideal monatomic gas is obtained as 



(i) 


where 7 = 5/3 is the adiabatic index, ks is the Boltzmann 
constant, ji — 0.61 is the mean particle weight and m p is 
the proton mass. We also show the characteristic mean free 
the gas, which is numerically defined as (e.g. 


A mfp -23kpc( io8 e K ) (10 -3 cm -3 ) ' (2) 

Notice, that the mean free path is below kpc in central few 
tens of kpc and increases with distance, reaching ~ 10 kpc 
at distance 400 kpc from the center. 


path A m f p in 


Sarazin[l988 
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Figure 2. Left: Chandra mosaic image of the Perseus Cluster in the 0.5 — 3.5 keV band. The units are counts/s/pixel. The size of each 
pixel is 1 arcsec. Right: residual image of the cluster (the initial image divided by the best-fitting spherically-symmetric / 3 — model of 
the surface brightness), which emphasizes the surface brightness fluctuations present in the cluster. Point sources are excised from the 
image. Black circles show the set of annuli used in the analysis of the fluctuations. The width of each annulus is 1.5' (~ 30 kpc). The 
outermost boundary is at a distance 10.5' 220 kpc) from the center. For display purposes, both images are lightly smoothed with a 

3" Gaussian. 



0 0.0007 0.0021 0.0049 0.01 0.022 0.044 0.088 0.18 0.35 0.7 )0e+00 2.89e+05 5.80e+05 8.70e+05 1.16e+06 1.45e+ 


Figure 3. Left: simulated map of the combined Chandra PSF within the mosaic image of the Perseus Cluster. The random positions of 
individual PSFs are used (see Section[2]for details). Right: the combined exposure map in seconds (lightly smoothed with a 3" Gaussian 
for display purposes). 
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Radius, kpc Radius, kpc 

Figure 4. Radial profiles of thermodynamic properties of the 
Perseus Cluster obtained from Chandra observations in the 0.6 — 9 
keV band. Orange (filled) and red (open) points with lcr error 
bars on the left panels show the deprojected gas number electron 
density n e , electron temperature T e and abundance of heavy el¬ 
ements relative to solar Z/Z§, assuming the latter one is con¬ 
stant 0.5 relative to Solar or is a free parameter in the fitting 
spectral model of an optically thin plasma respectively. Analytic 
approximations of the observed profiles by smooth functions are 
shown with black solid curves. Blue points on the top right panel 
show the spherically-symmetric X-ray surface brightness profile 
obtained from the 0.5 — 3.5 keV image, while the black curve 
shows its best-fitting (3 —model with the core radius r c = 1.26 
arcmin and the slope (3 = 0.53. The cluster sound speed c s and 
the mean free path Acalculated from the approximations are 
shown in the middle and bottom panels on the right. 


4 POWER SPECTRUM OF THE SURFACE 
BRIGHTNESS FLUCTUATIONS 

X-ray SB of galaxy clusters is typically peaked towards cen¬ 
ters and decreases rapidly with radius. This is especially 
true for relatively dynamically relaxed clusters like Perseus 
( Mantz et al.|2014 ). Therefore, the PS of the initial SB will 
be dominated by this gradient, leading to a large power on 
small wavenumbers k and contamination of the power on 
larger k. In order to avoid this, we first remove the global 
SB gradient and analyze fluctuating part of the SB only. 
In order to work with dimensionless units of the SB fluctua¬ 
tions, we divide the image of the Perseus Cluster by a simple 
spherically-symmetric [3 —model of the SE0 Fig .[2] shows the 
residual image of the SB fluctuations in Perseus. Notice a va- 


3 One can also subtract the model, however in this case addi¬ 
tional weights proportional to the global SB profile have to be 
introduced in order to treat the fluctuations at different radii 
equally. Clearly, the choice of this underlying model is arguable, 
but we will return to this problem in Section |6l 


riety of structures of different sizes and morphologies present 
in the cluster, especially within the central 3 arcmin. Let us 
quantify these structures statistically, by measuring their 
power as a function of length scale. 


We use a modified A—variance method (Arevalo et al. 


2012) to calculate the PS of SB fluctuations. This method is 


tuned to cope with non-periodic data with gaps, the PS of 
which are smooth functions. Dashed curves in Fig. [5] show 
the PS of SB fluctuations in the central 3 arcmin (central 
AGN is excluded) and in a broad annulus 3 — 10.5 arcmin. 
The spectra decrease with wavenumber k and reach constant 
value on small scales, where the signal is dominated by the 
Poisson noise. 

Knowing the number of counts n cn ts in each pixel in 
the cluster image in counts, the contribution of the Poisson 
noise is evaluated. A white noise has a flat spectrum; how¬ 
ever, small deviations are still possible due to imperfections 
of the A—variance method (e.g. due to filter behavior close 
to the edges of the image or near gaps, see | Arevalo et al.| 
(2012) for details). We made 100 realizations of the Pois¬ 
son noise simply multiplying y/n cn ts by a random number 
from a normal distribution (with mean 0 and variance 1) in 
each pixel. For each realization, the PS of the Poisson image 
is calculated and, finally, the mean and the scatter of the 
resulting spectrum are obtained. Dotted horizontal regions 
in Fig. [5] show the level of the Poisson noise with the un¬ 
certainties in both considered regions. Subtracting the PS 
of the noise from the total (SB fluctuations+noise) PS, the 
PS of SB fluctuations is obtained (hatched region with solid 
boundaries in Fig. [5). The uncertainties reflect the photon 
counting noise at high frequencies (lcr uncertainty is shown) 
and stochasticity of the signal at low frequencies. We do not 
separate these two types of uncertainties. Multiple realiza¬ 
tions of the white noise PS P wn estimate the scatter cr Pwn . 
These errors are then scaled to the measured PS of SB fluc¬ 
tuations P by a simple multiplication 


a p — cr Pwn . (3) 

-Mvn 

Based on experiments with the whole process pipeline, we 
conservatively added 3 per cent systematic uncertainty to 
the scatter of the Poisson noise PS, which takes into account 
subtle numerical and instrumental effects. 

Finally, the PS of SB fluctuations have to be corrected 
for the Chandra PSF and for the contribution of the un¬ 
resolved point sources on small scales. Performing multiple 
realizations of the PSF maps (see Section |2j), and each time 
calculating the PS of the map within the region of our inter¬ 
est, the mean spectrum of the PSF is obtained. Dash-dotted 
curves in Fig. [5] show the PSF PS in central 3 arcmin and 
in the annulus 3 — 10.5 arcmin. The deviations of the PSF 
shape from a S —function are reflected in a drop of power 
on the high-/c end of the spectrum. Dividing the PS of the 
SB fluctuations by the PSF power spectrum, the SB power 
suppressed by the PSF is recovered. 

The contribution of unresolved point sources is esti¬ 
mated by using the Perseus sensitivity map and by obtaining 
the shape and normalization of the Log N-Log S distribution 
of the resolved sources. In case of the Perseus Cluster, the 
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Figure 5. Left: power spectra (PS) of the X-ray surface brightness (SB) fluctuations in the Perseus Cluster in the central 3 arcmin 
(purple) and in the annulus 3 — 10.5 arcmin (orange). Dash curves: the initial PS. Dotted regions: PS of the Poisson noise. Dash-dotted 
curves: PS of the combined Chandra PSF in both analyzed regions. Hatched region with solid boundaries: PS of the SB fluctuations 
after subtraction of the Poisson noise. Both statistical (lcr) and stochastic uncertainties are taken into account. Notice, that very high 
statistics of counts accumulated by the 1.4 Ms observations allows us to probe fluctuations on a broad range of length scales, from > 100 
kpc down to few kpc (comparable with the mean free path). Right: amplitude of the SB fluctuations in the same regions, obtained from 
the PS corrected for the PSF (see Section [4|. Notice, that the amplitude of fluctuations varies with the radius. For these two regions a 
maximal difference by a factor of 4 is reached on ~ 17 kpc scale. 

contribution of the unresolved sources is negligible (within 
the cluster core). 

After applying all corrections described above, the char¬ 
acteristic amplitude A^d (k) of the SB fluctuations is calcu¬ 
lated, namely A 2 £>(/c) = ^/P2D(k)27vk 2 , where P 2 £>(/c) is the 
measured PS. The amplitude is a more convenient charac¬ 
teristic since its units are the same of the variable in a real 
space. Fig. [5] shows the amplitude of SB fluctuations in the 
central 3 arcmin region and in the outer annulus 3 — 10.5 
arcmin. High statistics of counts allows us to probe the SB 
fluctuations on a broad range of length scales - more than 
an order of magnitude, down to the scales comparable with 
the mean free path (see Fig. |4j). The amplitude varies from 
~ 3 — 14 per cent on scales 2 — 50 kpc in the central 3 arcmin 
to - 2 - 33 per cent on scales ~ 14 — 170 kpc in the 3 — 10.5 
arcmin region. The fact that the amplitudes in both regions 
differ at least by a factor of 3 over a broad range of scales 
shows that the amplitude of SB fluctuations in Perseus varies 
with radius. Therefore, instead of global characterization of 
the fluctuations in the cluster core, we consider below the 
fluctuations in a set of narrow radial annuli. 


5 RESULTS 

5.1 Power spectrum of density fluctuations 

We divided the cluster into a set of radial annuli with a 
width 1.5' (w 31 kpc) (Fig. [ 2 ]) and calculated the PS of SB 


fluctuations in each annulus, subtracting the Poisson noise 
and correcting for the unresolved point sources and the PSF. 
The latter one was done using PSF PS shown in Fig. [6] 
The fact that the X-ray SB lx ~ f n^dl for hot clusters 
with T e > 2 keV (see Fig. [TJ allows us to convert the 2D 
PS of SB fluctuations into the 3D PS of density fluctuations 
on scales smaller than the length of the line of sight. The 


procedure is described in Churazov et al. (2012). Namely, 


knowing the global spherically-symmetric model of the gas 
emissivity, the conversion factor between the 3D and 2D 
power spectra at each line of sight z is 


P2p(k) 

PsD(k) 


[ I\W(k z )\‘ 


dk z 


(4) 


where W(k z ) is the PS of the normalized emissivity distri¬ 
bution along the line of sight. 

We tested this procedure using high-resolution cosmo¬ 


logical simulations of galaxy clusters (Nagai, Vikhlinin, & 
Kravtsov|2007 Nelson et al.|2014 ). Using non-radiative sim¬ 
ulations, we obtained the PS of density fluctuations directly 
from the simulated 3D data for a sample of 6 relatively re¬ 
laxed clusters and recovered the PS from projected images 
of the X-ray SB (applying the same procedure we use for the 
observed X-ray images). The details of the sample and the 
analysis are described in |Zhuravleva et al. (2014a). Direct 
comparison of both PS (or amplitudes) confirms that the X- 
ray images can provide reasonably accurate measurements 


k, kpc -1 

10" 2 10- 1 1 
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Figure 6. Power spectrum (PS) of the combined Chandra PSF 
obtained in each annulus of the Perseus Cluster in our analysis 
(Fig# The combined PSF accounts for the shape of the PSF at 
different offsets and exposure of each observation (see Section [2J. 
Each spectrum is averaged over 10 spectra of different realizations 
of the PSF map. Since the PSF is not a 5—function, the PS deviate 
significantly from unity (dotted line), especially on scales smaller 
than 15 kpc (k 0.02 arcsec -1 ). 



Figure 7. Amplitude of density fluctuations in the inner 200 kpc 
(radius) region in the simulated relaxed galaxy cluster (CL21 clus¬ 
ter with a total mass M500c = 6.08 • 1O 14 M0 at rsooc = 1215.2 
kpc, see Zhuravleva et al.| ( |2014a| > for details) obtained from the 
3D density information (dash navyblue curve) and recovered from 
the image of the X-ray surface brightness (hatched orange/red re¬ 
gion). The width of the hatched region reflects the uncertainty of 
the recovered 3D amplitude of density fluctuations due to vari¬ 
ations of the converting geometrical factor with the projected 
radius. Notice that the true (from the 3D data) amplitude of 
density fluctuations is within the hatched region. This confirms 
that the observational procedure of measuring the amplitude of 
density fluctuations from the X-ray images recovers reasonably 
accurate the amplitude of the 3D density fluctuations. 


of the PS of density fluctuations. Fig. [ 7 ] shows an example 
of one of the clusters in our sample. 

The accuracy of the recovery of the 3D spectrum in 
each annulus depends on how steep the emissivity profile 
is. In case of the Perseus Cluster, the emissivity is strongly 
peaked towards the center, therefore the conversion factor 
significantly varies within the radial annuli. The narrower 
the annulus, the smaller the uncertainty. To do the conver¬ 
sion, we use the mean (within the annulus) value of this 
factor. In Section [6] we will discuss the uncertainties associ¬ 
ated with this step of the analysis. 

Fig. [5] shows the amplitude of density fluctuations A$d 
obtained in each annulus from the PS of SB fluctuations P 2 D 
as a function of wavenumber k. Scales, which are smaller 
than the width of each annulus (« 31 kpc) are plotted. 
Hatched regions show the amplitude on wavenumbers, over 
which we deem the measurements are least affected by the 
uncertainties and more robust against the observational lim¬ 
itations. The choice of the high-fc limit is set by the statisti¬ 
cal uncertainty or by the PSF distortions of the amplitude. 
For both cases, the uncertainty is less than 20 per cent in 
this hatched region. At low k, we exclude scales, on which 
the amplitude flattens or decreases towards smaller k. This 
flattening is most likely determined by the presence of sev¬ 
eral characteristic length scales (e.g. the distance from the 
center, scale heights, additional driver of fluctuations) and 
by the uncertainties in the choice of the underlying cluster 
model (see Section [6|. 

The amplitude of density fluctuations is less than 15 per 
cent, except for the central ~ 30 kpc. Table ^ summarizes 
the characteristic values of the amplitude in all considered 
regions. The amplitude in the innermost region is noticeably 
higher and the curve itself is flatter than in other regions. 
This is not surprising, since a significant fraction of the vol¬ 
ume of this region is dominated by sharp edges around bub¬ 
bles of relativistic plasma, shocks around them, filamentary 
structures and absorption features, flattening the high-A; end 
of the amplitude. Simple tests showed the drop of power 
on small scales when these features are excluded from the 
analysis. The radial variations of the amplitude at different 
scales is shown in the right panel of Fig. [8] Notice, that the 
amplitude decreases with the distance from the center and 
it increases with the scale (characteristic size) of the fluctu¬ 
ations. 


5.2 Velocity power spectrum 

Knowing the amplitude of density fluctuations on each 
length scale, the amplitude (and the PS) of velocities of gas 
motions present in the ICM can be obtained. Simple theo¬ 
retical arguments show that in the stratified atmospheres of 
relaxed galaxy clusters, the amplitude of density fluctuations 

A 3 o(k) = and one-component velocity V\ k are propor- 

po 

tional to each other at each length scale l = 1/k within the 
intertial range of scales, namely 
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Figure 8. Left: amplitude of density fluctuations in the Perseus Cluster versus wavenumber calculated in a set of radial annuli shown in 
Fig. [2] Hatched regions show the amplitude on scales where the measurements are least affected by systematic and statistical uncertainties. 
The width of each curve reflects estimated lcr statistical and stochastic uncertainties. Right: radial profiles of the amplitude of density 
fluctuations measured on certain scales (see the legend). Notice, that the amplitude decreases with the radius and increases with the 
length scale of fluctuations. 
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Figure 9. Left: amplitude of one-component velocity of gas motions versus wavenumber k = l/l , measured in a set of radial annuli (see 
the legend) in the Perseus Cluster. The velocity at each scale is obtained from the amplitude of density fluctuations, shown in Fig. [8] 
using relatio n [5] The color-coding and notations are the same as in Fig. [8] The slope of the amplitude for pure Kolmogorov turbulence 
(Kolmogorov 1941), k ~ 1 / 3 , is shown with dash line. Right: radial profiles of one-component velocity amplitude measured on certain 
length scales written in the legend. 
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Table 1. The mean values of the amplitude of density fluctuations and the amplitude of one-component velocity in the Perseus Cluster 
(see Fig. [5] and |9j . The values are given on scales l = 1/k: 10, 20 and 30 kpc. 


Annulus 

Amplitude of density fluctuations 
A 3 D = Sp/p 

Amplitude of 

one-component velocity 
Vi,fc, km/s 


l = 10 kpc 

1 = 20 kpc 

1 = 30 kpc 

l = 10 kpc 

l = 20 kpc 

1 = 30 kpc 

O' - 1.5' (0 - 31 kpc) 

0.22 

- 

- 

200 

- 

- 

1.5'- 3'(31- 62 kpc) 

0.11 

0.14 

- 

110 

140 

- 

3' - 4.5' (62 - 94 kpc) 

0.08 

- 

- 

90 

- 

- 

4.5'- 6'(94- 125 kpc) 

0.09 

0.1 

0.11 

110 

130 

140 

6 ' - 7.5' (125 - 156 kpc) 

- 

0.09 

0.12 

- 

120 

150 

7.5'- 9' (156 - 187 kpc) 

- 

0.08 

0.09 

- 

100 

120 

9'- 10.5' (187 - 219 kpc) 

- 

- 

0.08 

- 

- 

100 



Figure 10. Characteristic length scales present in the Perseus 
Cluster at different distances R from the center. Colored shaded 
regions: range of scales, on which our measurements are least 
affected by systematic and statistical uncertainties (hatched re¬ 
gions in Fig. [ 8 ] and Fig.§. Solid curve: the Ozmidov scale lo of 
turbulence in the core of Perseus, obtained assuming local bal¬ 
ance between turbulent heating and radiative cooling (relation 
[ 6 j . We do not plot lo in central ~ 20 kpc, since the measured gas 
entropy is flat there, leading to lo —> oo. Dash curve: the Kol¬ 
mogorov (dissipation) scale Ik for unmagnetized plasma (relation 
|7| . Dot-dash curve: the electron mean free path for unmagnetized 
plasma. See Section [5^2] for discussion. 


where po is the mean gas density and 771 is the proportion¬ 
ality coefficient (Zhuravleva et al. 2014a). If the injection 
scale of the turbulence is larger than or comparable with 
the Ozmidov scale ( Ozmidov|19 92; Breth ouwer et al.|2007 ), 
the coefficient 771 is set by gravity-wave physics on large, 
buoyancy-dominated scales, and is ~ 1 in atmospheres of 
galaxy clusters. It remains the same on smaller scales (within 
the inertial range) where the density becomes a passive 
scalar. Cosmological simulations of galaxy clusters confirm 


the relation |5] giving the averaged over a sample of relaxed 
clusters value 771 = 1 zb 0.3 ( Zhuravleva et al.|2014a ). 


Fig. [9] shows the amplitude of one-component veloc¬ 
ity versus wavenumber obtained from the relation [5] in all 
seven radial annuli in the Perseus core, as well as the ra¬ 
dial profiles of the amplitude on certain scales. Like in Fig. 
[ 8 ] hatched regions show the amplitude on scales, where the 
measurements are least affected by systematic and statistical 
uncertainties. Notice that the velocity amplitude is higher 
towards the center, suggesting a power injection from the 
center. The spectra (amplitudes) show similar dependence 
on k, larger fluctuations are on larger k , consistent with 
cascade turbulence. The slope of the velocity PS is broadly 
consistent with the slope for canonical Kolmogorov turbu¬ 
lence, accounting for the errors and uncertainties. Within 
the “robust” range of scales, Vi,k varies from ~ 70 km/s up 
to - 210 km/s on scales ~ 5 — 30 kpc (see Table [I]). The ve¬ 
locity amplitudes quantitatively match our expectations of 
typical velocities in the ICM from various observational con¬ 
straints (see e.g. |Churazov et al.|2004||Schuecker et al.|2004| 


Werner et al.||20()9| |Sanders et al.||20ib| |Sanders, Fabian, fc 

Smith |2011||de Plaa et al.||2012| |Sanders et al.||2013[ IPinto' 


tions (see e.g. Norman V Bryan 

1999 Dolag et al.| 

2005 

Iapichino V Niemeyer||2008| jLau, Kravtsov, V Nagai 

2009 

Vazza et al. 2011; Miniati|2014| and references therein). Fu- 


ture direct measurements of the velocity of gas motions with 
X-ray calorimeter on-board Astro-if observatory (Takahashi 


et al.|2014 ) will allow us to test the density-velocity relation 


5] Recently-proposed observational strategy for the Perseus 
Cluster maps the cluster core in radial and azimuthal di¬ 


rections, enabling us to do the calibration (Kitayama et al 


2014). 


It was recently shown that in the cores of Perseus and 
Virgo clusters, where the cooling time is shorter than the 
Hubble time, the heating of the gas due to dissipation of 
turbulence is sufficient to offset radiative cooling losses (Zhu- 
ravleva et al.||2014b ). Assuming that the dissipation rate is 


matching the cooling rate, it is straightforward to estimate 
the Ozmidov scale lo of the turbulence using only thermo¬ 
dynamic properties of the Perseus Cluster. Namely, 
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Figure 11. Left: residual image of the Perseus Cluster in the 0.5 —3.5 keV band (the same as Fig. [ 2 ]) with excluded bubbles of relativistic 
plasma taken from Fabian et al. ([2011) (black ovals). Right: one-component velocity amplitude versus wavenumber obtained from the 
image on the left. The color-coding and labels are the same as in Fig. [9] 


lo 


_ /v - 3/2 1/2 _ 


-3/2 ( Qcool \ 

V P ) 


1/2 


( 6 ) 


where N — 


lH s 


is the Brunt-Vaisala frequency in the 


cluster atmosphere, g is the acceleration of gravity, H s = 

( ] is the entropy scale height, and e is the dissipa- 

V dr J 

tion rate per unit mass. Here we assume pe ~ the cooling 
rate Q CO oi = n e riiA n (T ), where n e and rii are the number 
densities of electrons and ions, respectively, and A n (T) is 
the normalized gas cooling function (Sutherland V Dopita 


1993). 


Fig. [TO] shows the radial profile of the Ozmidov scale lo 
and a range of scales we are probing in each annulus (hatched 
regions in Fig. [ 9 ]). One can see that lo is within the interval 
of scales we are probing at each distance from the center 
within the cluster core. This means that the main assump¬ 


tion behind the derivation in Zhuravleva et al. (2014a) is 


satisfied, namely the Ozmidov scale is smaller than the in¬ 
jection scale of turbulence (assuming that we are probing 
velocity PS within the inertial range). Notice, that we do 
not show l 0 at R < 20 kpc since the measured gas entropy 
is fiat towards the center, leading to lo 00 . 

It is also interesting to compare the scales we are prob¬ 
ing with the Kolmogorov (dissipation) scale 


3/4 


Ik — 


{Qcool/p) 1 '^ 


where Vkin — 


L'dyn 


(7) 


is the kinematic viscosity, which is 


obtained through the dynamic viscosity Vdyn for an ion¬ 
ized plasma without magnetic field. Fig. [To] shows the Kol¬ 
mogorov scale Ik as well as the mean free path for compar¬ 


ison. The Kolmogorov scale is significantly below the scales 
we are probing, which justifies even better our assumption 
about the inertial range of scales. In reality, the situation 
can be more complicated, e.g., due to the presence of mag¬ 
netic fields. Here, we neglect these complications, following 
the simplest approach as a good starting point. 

We cannot claim, of course, that the structures in the 
SB are due to turbulence only. For example, there are many 
bubbles of relativistic plasma in the Perseus core (see e.g. 


Boehring er et al.||1993| |Churazov et al.||2000| |Fabian et al.| 


2000), which may contribute to the signal. The question is 
whether their contribution is dominant in the considered re¬ 
gions. Various sharp structures (e.g. bubble edges) would 
give the slope of the amplitude &; -1 / 2 , which is hard to dis¬ 
criminate from the Kolmogorov slope accounting for 

the uncertainties of our measurements and the assumptions 
used. However the contribution of the bubbles to the signal 
can be easily seen if we repeat the analysis excluding the 
known bubbles from the image of Perseus. Fig. |11| shows the 
one-component velocity amplitude calculated from such im¬ 
age. Notice, that the velocity amplitude decreases by a fac¬ 
tor ~ 1.6 — 1.2 depending on scale in the central 1.5 arcmin. 
This is expected since bubbles in this region are particularly 
prominent and occupy a substantial fraction of the volume. 
In the 1.5 — 3 arcmin annulus, this factor is ~ 1.15 over 
the whole range of scales. Outside 3 arcmin, the exclusion 
of bubbles from the analysis does not change the measured 
velocity amplitude. 


2009), mergers and gas sloshing (Markevitch V Vikhlinin 


Sound waves (Fabian et al. 2006 Sternberg V Soker 


2007) might also contribute to the observed spectrum of 
the density fluctuations. Unsharp masking of the Perseus 
Cluster revealed quasi-spherical structures (“ripples”) in the 
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SB, which have been interpreted as isothermal sound waves 
(Fabian et al. 2000 2003). An alternative interpretation 
is stratified turbulence, which arises naturally in the clus¬ 
ter atmosphere, where rough estimates give Froude number 
Fr ~ 0.3 — 1 (Zh uravleva et ah]|2014a ). In this case, the ra¬ 
dial size of each “ripple” A r is determined by HV/c s , where 
H is the characteristic scale height and V is the velocity 
amplitude. For example, in the 1.5 — 3 arcmin annulus in 
Perseus, the scale height H is ~ 40 — 70 kpc, the sound 
speed c s ~ 900 km/s and the velocity measured from the 
SB fluctuations is ~ 100 — 140 km/s. This gives us typical 
radial sizes of the fluctuations ~ 5 — 10 kpc, which are con¬ 
sistent with those seen in the X-ray image of the cluster once 
the large-scale asymmetry in SB is removed. Certainly, we 
cannot claim that the observed fluctuations are associated 
with the stratified turbulence only. The question on the na¬ 
ture of the fluctuations will be addressed in our future work. 
At the very least, the measured turbulent velocities can be 
treated as an upper limit. 


5.3 Gas clumping 


The term “gas clumping” refers to any deviations of gas 
density isosurfaces from equipotential surfaces and, there¬ 
fore, includes both large-scale inhomogeneities in the gas 
and small-scale clump^j It is difficult to determine or even 
define unambiguously the surfaces of constant potential in 
clusters. Therefore, in practice, we measure clumping rela¬ 
tive to a model, which describes a global gas distribution in 
the cluster potential well. Namely, the gas clumping factor 
is usually defined as 


r ((p/po) 2 ) 

( p/po > 2 


( 8 ) 


where () denotes the mean or median inside a spherical shell 
and po is the global density model, which accounts for the 
density gradient within each shell. 

Gas clumping leads to an overestimate of the gas density 
and, as a consequence, affects the gas mass and total hydro¬ 


static mass measurements of galaxy clusters (see e.g. Lau, 
Kravtsov, & Nagai 2009]) as well as the SZ measurements of 
the Compton parameter (e.g. Khedekar et al.|2 013). Numer¬ 
ical simulations predict clumping < 1.05 in central 0.5t*5oo, 


~ 1.1 — 1.4 at 7*500 and up to 2 at 7*200 (see e.g. 

Math- 

iesen, Evrard, & Mohr 1999; Nagai & Lau 2011 

Roncarelli 

et al.|2013| Zhuravleva et al. 2013a; |Vazza et al. 

2013). The 


clumping factor obtained from X-ray observations of cluster 
outskirts is slightly higher and varies from ~ 2 to ~ 3 at 
7*200 (see e.g. Urban et al.|2014 Walker et al.|2012 |Mo randi 


&; Cui|2014[ and references therein). The main challenges of 


the clumping measurements from the X-ray observations are 
a very low X-ray SB in cluster outskirts, where clumping is 
expected to be higher, and, in contrast, a low value of the 
clumping factor, less than 5 per cent according to numerical 
simulations, in the bright cluster cores. 


4 Less general definitions are often used in astrophysics. For ex¬ 
ample, clumping can be referred to gas clumps and subhalos only, 
which are notably denser and cooler than the ambient gas. 



Figure 12. Radial profile of the clumping factor C — 1 in the 
core of the Perseus Cluster obtained from the relation § us¬ 
ing the measurements of density fluctuations power spectra. Red 
points: the measured density spectrum is fitted by a power-law 
function with the Kolmogorov slope —11/3 and extended to small 
and large scales, which we do not probe in our analysis. Blue 
points: the slope of the power-law function is assumed —3.1 and 
—4 for the innermost annulus 0 — 1.5 and for other six regions re¬ 
spectively. The slope —4 is obtained from the power spectrum of 
density fluctuations in the broad annulus 1.5 — 10.5 arcmin. Dash 
points: k m i n = 1/R, kmax = oo. Solid points: k m i n = 0.01 
arcmin -1 (inverse width of each annulus), kmax = oo. Gray dot¬ 
ted points: clumping factor obtained from the measured power 
spectrum of density fluctuation directly, avoiding any approxi¬ 
mations with the power-laws. Measurements within the range 
of scales least affected by the uncertainties (hatched regions in 
Fig.§ are used. The corresponding scales are written below each 
point. Notice, that the clumping factor is less than 7 — 8 per cent 
outside the central ~ 30 kpc. 


Presenting the density p through unperturbed po and 
fluctuating Sp/po parts, p — po(l + ^p/po), the gas clumping 
factor definition ([8]) can be re-written through the PS of 
density fluctuations as 


C = 1 + 



1 + 


K ma. 

/ 


4tt P 3D (k)k 2 dk, 


(9) 


where the amplitude of density fluctuations — is defined 
2 P ° 

through the PS as ( — J = 


po 


J 4:7rP3D(k)k 2 dk and the in¬ 


tegration is over all wavenumbers, i.e. from kmin = 0 till 
kmax — oo. Using this relation and our measurements of 
the PS of density fluctuations in the Perseus core, we es¬ 
timate the gas clumping factor, which is shown in Fig. [12] 
For each annulus we obtain kmin and kmax as the smallest 
and the largest wavenumbers within the range of wavenum¬ 
bers least affected by uncertainties (hatched regions in Fig. 
|8|, see dotted gray points in Fig. 12 The gas clumping is 
very low, less than 2 per cent at R 30 kpc, which is 
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Figure 13. Top: underlying models of “unperturbed” surface brightness in the Perseus Cluster. From left to right: spherically- 
symmetrical (3 —model (default choice), patched (3 —models with o = 80 arcsec (removes large-scale asymmetry), o = 30 arcsec and 
o = 10 arcsec (see Sect ion [6.1| for the details). Bottom: residual images of the SB fluctuations in Perseus obtained from the initial image 
divided by the underlying model (top panels). The smaller the o the smaller the structures included to the model and the less structures 
remain in the residual image. 


expected since we are integrating over a narrow range of 
wavenumbers. In order to estimate the total clumping, we 
approximated the PS of density fluctuations with power-law 
functions and extended them up to kmax = oo and down to 
kmin — 0.01 arcmin -1 (inverse width of each annulus, solid 
points in Fig. 12) or down to kmin = 1/R (dash points in 
Fig. 12), where R is the distance from the center to each 
annulus. Red points in Fig. [12] show the gas clumping in 
case when we fit a power-law with the Kolmogorov slope 
— 11/3 to the measured spectra, varying only its normaliza¬ 
tion in each annulus. For the blue points, we assume the 
slope being consistent with the slope of the global density 
spectrum obtained in the broad 1.5 — 10.5 arcmin annulus. 
Namely, for annuli at R > 30 kpc (1.5 arcmin), the slope 
of the spectrum is —4, while for the innermost annulus, the 
slope is —3.1. Notice, that the gas clumping is dominated 
by the large scales; it decreases with the distance from the 
center if we do not account for fluctuations on scales larger 
than the width of each annulus (it simply reflects the fact 
that the amplitude of fluctuations decreases with the radius, 
see Fig. |8|; outside the central ~ 30 kpc, the gas clumping 
is lower than 7 — 8 per cent, which leads to a density bias 
~ (C — l)/2 less than 3 — 4 per cent. 


6 UNCERTAINTIES IN THE ANALYSIS 

In this Section, we examine various systematic uncertainties 
in the analysis. Namely, we consider uncertainties associated 
with 

• the choice of the unperturbed cluster density model; 

• the bias of the A—variance method used for the PS 
calculations; 

• inhomogeneous exposure coverage; 

• conversion P 2 d to P^d] 

• velocity measurements. 

Even though the list of uncertainties is, admittedly, long, a 
conservative estimate of the total uncertainty on the ampli¬ 
tude of density fluctuations (one-component velocity) is less 
than 50 (60) per cent. 


6.1 Underlying model of the surface brightness 

Decomposition of the cluster image into “perturbed” and 
“unperturbed” components is ambiguous. Physically mo¬ 
tivated, the underlying model of the “unperturbed” com¬ 
ponent should reflect the global potential of the cluster in 
equilibrium, i.e. when the isosurfaces of the gas density and 
temperature are aligned with the equipotential surfaces. Our 
default choice of the model - /3 —model for the azimuthally 
averaged SB (or just mean SB) - seems to be a reasonable 
choice of such “unperturbed” model at least for (or close to) 
relaxed galaxy clusters. However, one can see the large-scale, 
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Figure 14. Amplitude of density fluctuations in the 3—4.5 arcmin 
annulus in the Perseus Cluster obtained by using spherically- 
symmetrical (3 —model as the underlying one (red) and more flex¬ 
ible models shown in Fig. |13| Notice, that the removal of the 
global cluster asymmetry (blue hatched region) does not affect 
the amplitude measurements relative to the default model (red) 
on scales we are probing. Also, notice that even if the amplitude 
suppression is present on scales, which are included to the un¬ 
derlying model (as expected), the amplitude on smaller scales, 
which are not affected by the model, remains almost the same. 
This means that the measured amplitude is indeed due to the 
presence of fluctuations of these scales and not due to the power 
leakage from the large scales. 


west - east asymmetry in the SB of the Perseus Cluster (Fig. 
[2] left), which might question our most simple and conser¬ 
vative choice. To check how the global cluster asymmetry 
affects the density amplitude measurements, we repeated 
the analysis, using two other /3—models of the SB averaged 
in 180° sectors. The best-fitting parameters of the /3 —model 
in the east (90° 4-270°) and the west (—90° 4-90°) sectors are 
r c = 1.76 and 0.89 arcmin, f3 m. 0.61 and 0.48, respectively. 
Choosing these underlying models, the maximal deviations 
of the density amplitude from the default case are on the 
largest scales in each annulus (typically ~ 30 — 20 kpc) and 
are less than 10 per cent in the central 6 arcmin, less than 
17 per cent in the 6 — 7.5 arcmin annulus and less than 30 
per cent at the distance 7.5 — 10.5 arcmin from the cluster 
center. 

Going beyond this simple spherically-symmetric model 
implies that we believe that the underlying cluster potential 
is more intricate. It is not clear to what degree of complexity 
of the model we should go. There is always a danger that 
some of the structures unrelated to the cluster gravitational 
potential are removed. We made a number of tests, introduc¬ 
ing more flexibility to the /3— model, by patching it on large 
scales to account for possible complexity of the potential. 
Namely, our patched model is defined as I prn = IpS a [lx / Ip ], 
where Ip is the f3 —model of azimuthally-averaged SB, S cr [*] 
denotes Gaussian smoothing with the smoothing window 


size a and lx is the clusters X-ray SB. Varying the size of 
the smoothing window, the model changes from most conser¬ 
vative symmetrical [3 —model (large a) to most complicated 
(and clearly implausible) one, which accounts for structures 
on all scales (small a) FI Examples of such models and the 
corresponding residual images of the SB fluctuations are 
shown in Fig. 13 Models with large a (e.g. ~ 80 arcsec) 
remove the large-scale (east-west) asymmetry of the clus¬ 
ter, while those with a smaller a absorb more features of 
the image on smaller scales. Fig. [14] shows the amplitude 
of density fluctuations in 3 — 4.5 arcmin annulus measured 
from the residual images in Fig. [13] Notice, that the re¬ 
moval of the large-scale asymmetry (blue hatched region) 
does not change the amplitude of density fluctuations rela¬ 
tive to spherically-symmetric /3 —model (red hatched region) 
over the range of scales probed in our analysis (< 30 kpc). 
As expected, the more complex models (with the smaller a), 
the broader range of scales, on which the amplitude is sup¬ 
pressed, and the stronger the suppression of the amplitude 
on the largest scales. 

One can estimate the scales, on which the amplitude 
is expected to be suppressed depending on the size of the 
window function a used in the patched underlying models. 
Assuming that the PS of density fluctuations is a power-law 
k~ a , the convolution of the image PS with the Mexican Hat 
filter will give (see relation A8 in Arevalo et al. 2012) 

4 


-Pmh(kr) OC 


M 


— 2{k/k r ) in 


d U k , 


( 10 ) 


where k r is the characteristic wavenumber. The convolu¬ 
tion of the PS of the image, with removed large-scale part 
smoothed with a Gaussian, is 


-Pmh,er {kr ) CK) 


/*" 


0 -2(k/k r ) 2 


— 2-7T 2 fe 2 cr 2 


) 2 cTfc.(ll) 


The ratio of both will give us an estimate of a wavenum¬ 
ber, on which the amplitude will be suppressed for any size 
a of the window function used for the underlying model of 
the SB. For the Kolmogorov PS, /c -11 ^ 3 , the suppression of 
the amplitude obtained using the patched f3 —model relative 
to the amplitude measured using the spherically-symmetric 
model is ~ 20 per cent on km 0.8/cr. This is roughly consis¬ 
tent with what we see in Fig. |14| For example, the underlying 
patched model with a = 20 arcsec gives a 20 per cent differ¬ 
ence in the amplitudes of density fluctuations on k c h ar ~ 0.04 
arcsec -1 . On k < k c h ar the “patched” amplitude is strongly 
suppressed, while on k > k c har the amplitudes remain al¬ 
most the same. This means that at k > k c h ar the measured 
amplitude of density fluctuations is due to the presence of 
fluctuations on these scales and not due to the leakage of 
power from the large scales. The net conclusion is that the 
amplitude of density fluctuations measured on scales < 30 
kpc is almost not affected by the choice of the underlying 
model, unless the cluster potential is very disturbed. 


5 An alternative way to remove large-scale structure is used in 
the analysis of SB fluctuations in the AWM7 cluster (Sanders V 
Fabian 2012). The cluster is modeled by fitting ellipses to contours 
of SB, spaced logarithmically. 
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6.2 Bias in the A— variance method 


The normalization of the PS obtained through the 
A—variance method may be biased slightly, depending on 
the slope of the PS (see Appendix B in Arevalo et al.|2012 ). 
The approximations of the initial PS of SB fluctuations with 
a power-law functions give the slopes ~ — 3-3.8 depend¬ 

ing on the distance from the cluster center. Therefore, the 
normalization of the measured amplitude of density fluctua¬ 
tions is on average overestimated by ~ 20 — 30 per cent and 
slightly lower, by ~ 10 per cent, in the innermost (0 — 1.5 
arcmin) and outermost (9—10.5 arcmin) annuli. 


6.3 Inhomogeneous exposure coverage 

The exposure map is not uniform and the brightness of the 
cluster itself varies across each annulus as seen in Fig.[2]and 
[3] Both arguments may bias the measured amplitude of den¬ 
sity fluctuations. By default, when calculating the amplitude 
by taking RMS of fluctuations present in filtered images, we 
use most uniform weighting scheme, w = 1, i.e. we treat 
all pixels in the image with the same weight. The measured 
amplitude is then more sensitive to fluctuations close to the 
outer edge of each annulus. The least uniform scheme, but 
at the same time most optimal for the reduction of the Pois¬ 
son noise requires weights to be w\ oc t exp Imod , where t eX p 
is the exposure map and Imod is the underlying model of the 
SB. In this case, those parts of the cluster that have higher 
number of counts would have larger weights. We also exper¬ 
imented with the weight W 2 oc t exp , which is more sensitive 
to the deepest-exposure parts of the image. We find <15 per 
cent (< 30 per cent) higher (lower) value of the amplitude 
of density fluctuations on scales < 20 kpc if weights w\ and 
W 2 are used. On larger scales, ~ 30 kpc, the amplitude can 
be < 60 per cent lower. As an example, Fig. |15| shows the 
amplitude in 1.5 — 3 arcmin annulus measured using three 
different weighting schemes. 


6.4 Conversion from 2D to 3D power spectrum 

We convert the two-dimensional PS of the SB fluctuations 
P 2 D into the three-dimensional PS of density fluctuations 
P 3 D using relation ©• The conversion factor f 2 D^ 3 D = 
P 2 D/P 3 D depends on a distribution of fluctuations along 
the line of sight and its value varies with projected distance 
from the cluster center. If radial profile of the SB is steep, the 
uncertainty can be large unless narrow annuli are considered. 
We use the value of f 2 D^ 3 D evaluated at the mean radius 
of each annulus. This factor is different in inner and outer 
edges of each annulus, leading to a maximal uncertainty on 
the measured amplitude < 17 per cent. Fig. [15] shows the 
amplitude in 1.5 — 3 arcmin annulus obtained, accounting 
for the variations of the conversion factor f 2 D^ 3 D within 
the annulus. 


6.5 Uncertainties in velocity measurements 

The proportionality between the amplitude of density fluc¬ 
tuations and velocity with the coefficient 77 ~ 1 is based on 


0.05 k * k P c ~‘ 



k, arcsec -1 

Figure 15. Amplitude of density fluctuations in 1.5 — 3 arcmin 
annulus in the Perseus Cluster measured varying the weight (top 
panel, see Section ( 


panel, see Section 
hatched region. 


6.3) and the conversion factor f2D^3D (bottom 


6.4). The default choice is shown with navyblue 


the simplest approach which is a good approximation (as 
confirmed by numerical simulations) that captures the key 
physics. The derivation neglects gas heating, cooling, heat 
fluxes, magnetic fields and assume that all perturbations 
are small. Magnetic fields on large scales should not signifi¬ 
cantly modify buoyancy physics since [3 — SirnkT/B 2 » 1. 
On small scales, the density fluctuations would behave as a 
passive scalar (Schekochihin et al. 2009). Therefore, we ex¬ 
pect the linear relation to hold even in MHD case, however 
the proportionality coefficient is likely to change in this case. 
The Astro-H X-ray observatory will allow us to verify and 
calibrate the relation. Any deviation from it will indicate 
the importance of the neglected physics. 

The intermittency of density fluctuations and velocity 
can be admittedly large in the Perseus Cluster. For exam¬ 
ple, analyzing the fluctuations in small patches within the 
3 — 4.5 arcmin annulus, we noticed that the amplitude of 
density fluctuations varies in a relatively broad range from 
3 to 10 per cent. In order to achieve statistical convergence, 
we perform our measurements in relatively wide annuli. The 
amplitude obtained in the twice broader annuli is consistent 
with the amplitudes measured in two individual annuli. 


7 CONCLUSIONS 

We performed detailed analysis of the X-ray SB and gas 
density fluctuations in a set of radial annuli within the core 
(central r ~ 220 kpc) of the Perseus Cluster, using deep 
Chandra observations. 

To summarize our findings: 

• The characteristic amplitude of the density fluctuations 
© 2015 RAS, MNRAS 000, [TUffi] 
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varies from 8 to 12 per cent on scales ~ 10 — 30 kpc within 
30 — 160 kpc annulus and from 9 to 7 per cent on scales 
~ 20 — 30 kpc in the outer annuli, 160 — 220 kpc. The am¬ 
plitude in the innermost 30 kpc is higher, up to 20 — 22 
per cent on scales ~ 5 — 15 kpc. The higher amplitude in 
this region reflects the presence of bubbles, shocks and sound 
waves around them, filaments and absorption feature, which 
occupy a large area of the considered region. The smallest 
scale we probe varies from ~ 5 kpc in the central 30 kpc to 
~ 25 kpc at distance 200 kpc from the center (while the 
mean free path for unmagnetized plasma is ~ 0.1 and 5 kpc, 
respectively). 

• Given stratification and gas entropy gradient in the at¬ 
mosphere of the cluster, we use linear relation between the 
amplitude of density fluctuations and velocity of gas motions 
to evaluate the characteristic velocity amplitude on different 
scales. The typical amplitude of the one-component velocity 
outside the central 30 kpc region is ~ 90 — 140 km/s on 
~ 20 — 30 kpc scale and ~ 70 — 100 km/s on smaller scales 
~ 7—10 kpc. These measurements match our expectations of 
typical velocities in the ICM from numerical simulations and 
various observational constraints. Measured velocity spectra 
suggest power injection from the center (e.g. from the cen¬ 
tral AGN in Perseus). Spectra are consistent with the cas¬ 
cade turbulence. Their slopes are broadly consistent with 
the slope for the canonical Kolmogorov turbulence. It was 
previously shown that the heating of the gas due to dissi¬ 
pation of such motions balances the gas radiative cooling in 
the cluster core ( Zhuravleva et ah||2014b ). 

• The gas clumping estimated from the PS of the density 
fluctuations is lower than 7 — 8 per cent at distance from 30 
to 220 kpc from the center, which gives a density bias less 
than 3 — 4 per cent in the cluster core. The clumping factor 
is dominated by fluctuations on large scales. 

• Systematic uncertainties in the analysis were analyzed. 
Conservative estimates of the final uncertainty on the am¬ 
plitude of density fluctuations is ~ 50 per cent and on the 
velocity amplitude is ~ 60 per cent. 

Future direct measurements of the velocities of gas mo¬ 
tions with the X-ray microcalorimeters on-board Astro-H , 
Athena and Smart-X will allow us to better calibrate the 
statistical relation between the amplitude of density fluc¬ 
tuations and the velocity amplitude. Any strong deviations 
from the proportionality coefficient ~ 1 would indicate the 
importance of the neglected physics in the ICM. 
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